An Artificial Intelligence Generated Automated Algorithm to Measure Total Kidney Volume in ADPKD.

Introduction: Accurate tools to inform individual prognosis in patients with autosomal dominant polycystic kidney disease (ADPKD) are lacking. Here, we report an artificial intelligence (AI)-generated method for routinely measuring total kidney volume (TKV). Methods: An ensemble U-net algorithm was created using the nnUNet approach. The training and internal cross-validation cohort consisted of all 1.5T magnetic resonance imaging (MRI) data acquired using 5 different MRI scanners (454 kidneys, 227 scans) in the CYSTic consortium, which was first manually segmented by a single human operator. As an independent validation cohort, we utilized 48 sequential clinical MRI scans with reference results of manual segmentation acquired by 6 individual analysts at a single center. The tool was then implemented for clinical use and its performance analyzed. Results: The training or internal validation cohort was younger (mean age 44.0 vs. 51.5 years) and the female-to-male ratio higher (1.2 vs. 0.94) compared to the clinical validation cohort. The majority of CYSTic patients had PKD1 mutations (79%) and typical disease (Mayo Imaging class 1, 86%). The median DICE score on the clinical validation data set between the algorithm and human analysts was 0.96 for left and right kidneys with a median TKV error of -1.8%. The time taken to manually segment kidneys in the CYSTic data set was 56 (±28) minutes, whereas manual corrections of the algorithm output took 8.5 (±9.2) minutes per scan. Conclusion: Our AI-based algorithm demonstrates performance comparable to manual segmentation. Its rapidity and precision in real-world clinical cases demonstrate its suitability for clinical application.

South Asian medical cohorts reveal strong founder effects and high rates of homozygosity.

The benefits of large-scale genetic studies for healthcare of the populations studied are well documented, but these genetic studies have traditionally ignored people from some parts of the world, such as South Asia. Here we describe whole genome sequence (WGS) data from 4806 individuals recruited from the healthcare delivery systems of Pakistan, India and Bangladesh, combined with WGS from 927 individuals from isolated South Asian populations. We characterize population structure in South Asia and describe a genotyping array (SARGAM) and imputation reference panel that are optimized for South Asian genomes. We find evidence for high rates of reproductive isolation, endogamy and consanguinity that vary across the subcontinent and that lead to levels of rare homozygotes that reach 100 times that seen in outbred populations. Founder effects increase the power to associate functional variants with disease processes and make South Asia a uniquely powerful place for population-scale genetic studies.

Frequent LPA KIV-2 Variants Lower Lipoprotein(a) Concentrations and Protect Against Coronary Artery Disease.

BACKGROUND: Lipoprotein(a) (Lp(a)) concentrations are a major independent risk factor for coronary artery disease (CAD) and are mainly determined by variation in LPA. Up to 70% of the LPA coding sequence is located in the hypervariable kringle IV type 2 (KIV-2) region. It is hardly accessible by conventional technologies, but may contain functional variants. OBJECTIVES: This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD. METHODS: We genotyped 4733G>A in the GCKD (German Chronic Kidney Disease) study (n = 4,673) by allele-specific polymerase chain reaction, performed minigene assays, identified proxy single nucleotide polymorphisms and used them to characterize its effect on CAD by survival analysis in UK Biobank (n = 440,234). Frequencies in ethnic groups were assessed in the 1000 Genomes Project. RESULTS: The 4733G>A variant (38.2% carrier frequency) was found in most isoform sizes. It reduces allelic expression without abolishing protein production, lowers Lp(a) by 13.6 mg/dL (95% CI: 12.5-14.7; P < 0.0001) and is the strongest variance-explaining factor after the smaller isoform. Splicing of minigenes was modified. Compound heterozygosity (4.6% of the population) for 4733G>A and 4925G>A, another KIV-2 splicing mutation, reduces Lp(a) by 31.8 mg/dL and most importantly narrows the interquartile range by 9-fold (from 42.1 to 4.6 mg/dL) when compared to the wild type. In UK Biobank 4733G>A alone and compound heterozygosity with 4925G>A reduced HR for CAD by 9% (95% CI: 7%-11%) and 12% (95% CI: 7%-16%) (both P < 0.001). Frequencies in ethnicities differ notably. CONCLUSIONS: Functional variants in the previously inaccessible LPA KIV-2 region cooperate in determining Lp(a) variance and CAD risk. Even a moderate but lifelong genetic Lp(a) reduction translates to a noticeable CAD risk reduction.

Discontinuation versus continuation of renin-angiotensin-system inhibitors in COVID-19 (ACEI-COVID): a prospective, parallel group, randomised, controlled, open-label trial.

BACKGROUND: SARS-CoV-2 entry in human cells depends on angiotensin-converting enzyme 2, which can be upregulated by inhibitors of the renin-angiotensin system (RAS). We aimed to test our hypothesis that discontinuation of chronic treatment with ACE-inhibitors (ACEIs) or angiotensin II receptor blockers (ARBs) mitigates the course o\f recent-onset COVID-19. METHODS: ACEI-COVID was a parallel group, randomised, controlled, open-label trial done at 35 centres in Austria and Germany. Patients aged 18 years and older were enrolled if they presented with recent symptomatic SARS-CoV-2 infection and were chronically treated with ACEIs or ARBs. Patients were randomly assigned 1:1 to discontinuation or continuation of RAS inhibition for 30 days. Primary outcome was the maximum sequential organ failure assessment (SOFA) score within 30 days, where death was scored with the maximum achievable SOFA score. Secondary endpoints were area under the death-adjusted SOFA score (AUCSOFA), mean SOFA score, admission to the intensive care unit, mechanical ventilation, and death. Analyses were done on a modified intention-to-treat basis. This trial is registered with ClinicalTrials.gov, NCT04353596. FINDINGS: Between April 20, 2020, and Jan 20, 2021, 204 patients (median age 75 years [IQR 66-80], 37% females) were randomly assigned to discontinue (n=104) or continue (n=100) RAS inhibition. Within 30 days, eight (8%) of 104 died in the discontinuation group and 12 (12%) of 100 patients died in the continuation group (p=0·42). There was no significant difference in the primary endpoint between the discontinuation and continuation group (median [IQR] maximum SOFA score 0·00 (0·00-2·00) vs 1·00 (0·00-3·00); p=0·12). Discontinuation was associated with a significantly lower AUCSOFA (0·00 [0·00-9·25] vs 3·50 [0·00-23·50]; p=0·040), mean SOFA score (0·00 [0·00-0·31] vs 0·12 [0·00-0·78]; p=0·040), and 30-day SOFA score (0·00 [10-90th percentile, 0·00-1·20] vs 0·00 [0·00-24·00]; p=0·023). At 30 days, 11 (11%) in the discontinuation group and 23 (23%) in the continuation group had signs of organ dysfunction (SOFA score ≥1) or were dead (p=0·017). There were no significant differences for mechanical ventilation (10 (10%) vs 8 (8%), p=0·87) and admission to intensive care unit (20 [19%] vs 18 [18%], p=0·96) between the discontinuation and continuation group. INTERPRETATION: Discontinuation of RAS-inhibition in COVID-19 had no significant effect on the maximum severity of COVID-19 but may lead to a faster and better recovery. The decision to continue or discontinue should be made on an individual basis, considering the risk profile, the indication for RAS inhibition, and the availability of alternative therapies and outpatient monitoring options. FUNDING: Austrian Science Fund and German Center for Cardiovascular Research.

Contamination detection in sequencing studies using the mitochondrial phylogeny.

Within-species contamination is a major issue in sequencing studies, especially for mitochondrial studies. Contamination can be detected by analyzing the nuclear genome or by inspecting polymorphic sites in the mitochondrial genome (mtDNA). Existing methods using the nuclear genome are computationally expensive, and no appropriate tool for detecting sample contamination in large-scale mtDNA data sets is available. Here we present haplocheck, a tool that requires only the mtDNA to detect contamination in both targeted mitochondrial and whole-genome sequencing studies. Our in silico simulations and amplicon mixture experiments indicate that haplocheck detects mtDNA contamination accurately and is independent of the phylogenetic distance within a sample mixture. By applying haplocheck to The 1000 Genomes Project Consortium data, we further evaluate the application of haplocheck as a fast proxy tool for nDNA-based contamination detection using the mtDNA and identify the mitochondrial copy number within a mixture as a critical component for the overall accuracy. The haplocheck tool is available both as a command-line tool and as a cloud web service producing interactive reports that facilitates the navigation through the phylogeny of contaminated samples.

Investigation of a nonsense mutation located in the complex KIV-2 copy number variation region of apolipoprotein(a) in 10,910 individuals.

BACKGROUND: The concentrations of the highly atherogenic lipoprotein(a) [Lp(a)] are mainly genetically determined by the LPA gene locus. However, up to 70% of the coding sequence is located in the complex so-called kringle IV type 2 (KIV-2) copy number variation, a region hardly accessible by common genotyping and sequencing technologies. Despite its size, little is known about genetic variants in this complex region. The R21X variant is a functional variant located in this region, but it has never been analyzed in large cohorts. METHODS: We typed R21X in 10,910 individuals from three European populations using a newly developed high-throughput allele-specific qPCR assay. R21X allelic location was determined by separating the LPA alleles using pulsed-field gel electrophoresis (PFGE) and typing them separately. Using GWAS data, we identified a proxy SNP located outside of the KIV-2. Linkage disequilibrium was determined both statistically and by long-range haplotyping using PFGE. Worldwide frequencies were determined by reanalyzing the sequencing data of the 1000 Genomes Project with a dedicated pipeline. RESULTS: R21X carriers (frequency 0.016-0.021) showed significantly lower mean Lp(a) concentrations (- 11.7 mg/dL [- 15.5; - 7.82], p = 3.39e-32). The variant is located mostly on medium-sized LPA alleles. In the 1000 Genome data, R21X mostly occurs in Europeans and South Asians, is absent in Africans, and shows varying frequencies in South American populations (0 to 0.022). Of note, the best proxy SNP was another LPA null mutation (rs41272114, D' = 0.958, R2 = 0.281). D' was very high in all 1000G populations (0.986-0.996), although rs41272114 frequency varies considerably (0-0.182). Co-localization of both null mutations on the same allele was confirmed by PFGE-based long-range haplotyping. CONCLUSIONS: We performed the largest epidemiological study on an LPA KIV-2 variant so far, showing that it is possible to assess LPA KIV-2 mutations on a large scale. Surprisingly, in all analyzed populations, R21X was located on the same haplotype as the splice mutation rs41272114, creating "double-null" LPA alleles. Despite being a nonsense variant, the R21X status does not provide additional information beyond the rs41272114 genotype. This has important implications for studies using LPA loss-of-function mutations as genetic instruments and emphasizes the complexity of LPA genetics.

Reference-based phasing using the Haplotype Reference Consortium panel.

Haplotype phasing is a fundamental problem in medical and population genetics. Phasing is generally performed via statistical phasing in a genotyped cohort, an approach that can yield high accuracy in very large cohorts but attains lower accuracy in smaller cohorts. Here we instead explore the paradigm of reference-based phasing. We introduce a new phasing algorithm, Eagle2, that attains high accuracy across a broad range of cohort sizes by efficiently leveraging information from large external reference panels (such as the Haplotype Reference Consortium; HRC) using a new data structure based on the positional Burrows-Wheeler transform. We demonstrate that Eagle2 attains a ∼20× speedup and ∼10% increase in accuracy compared to reference-based phasing using SHAPEIT2. On European-ancestry samples, Eagle2 with the HRC panel achieves >2× the accuracy of 1000 Genomes-based phasing. Eagle2 is open source and freely available for HRC-based phasing via the Sanger Imputation Service and the Michigan Imputation Server.